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We show that the addition of a magnetic field parallel to the current induces self sustained intrin- 
sic current oscillations in an asymmetric double barrier structure. The oscillations are attributed 
to the nonlinear dynamic coupling of the current to the charge trapped in the well, and the effect 
of the external field over the local density of states across the system. Our results show that the 
system bifurcates as the field is increased, and may transit to chaos at large enough fields. 
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Since the first observation of resonant tunneling 
through a semiconductor double barrier structure (DBS) 
by Chang, Esaki and Tsu, devices based on this struc- 
ture have been fou»d to exhibit a variety of interesting 
physical properties^. Besides a peak due to simplejeso- 
nant transmission- structures arising from phonona and 
plasmon -assistedu resonant tunneling, as well as from 
Landau level matchingoB, have been observed. More- 
over, an intrinsic dynamical bistability and hysteresis 
in the negative differential resistance region (NDRLof 
the current- voltage characteristic has been reportedo#Q. 
This is a non-linear effect produced by the Coulomb re- 
pulsion experienced by the incoming electrons fmm the 
charge buildup in the space between the barriers^. 

Qualitative arguments given by Ricco and AzbelH sug- 
gest the existence of intrinsic oscillations in a DBS for the 
one dimensional case. When the energy of the incoming 
electrons matches the resonance energy they enter the 
device and charge the potential well, lifting its bottom, 
thus driving the system away from resonance. The ensu- 
ing current decrease lowers the charge in the well bring- 
ing the system back to resonance, and a new cycle in an 
oscillatory behavior begins. Numerical calculations sup- 
port this prediction for ballistic electrons provided they 
are not monoenergetic, and in fact their energy [Spread is 
comparable, or wider than the resonance widthEI. 

In this letter we report the presence of sustained intrin- 
sic oscillations when the source is a Fermi sea of electrons, 
as one normally has in DBS devices. The oscillations 
may be regular or irregular, and are made possible by 
the presence of a magnetic field applied in the direction 
of the current. They arise thanks to the drastic change 
that the field causes on the density of states, and from the 
magnetic field induced enhancement of the nonlinearity 
of the system, as described below. 

Consider a DBS device. In order to study its time 
evolution under a bias we adopt a first-neighbors tight- 
binding model for the Hamiltonian. Inclusion of a mag- 
netic field B in the growth direction, henceforth called 
the z direction, is simple if a parabolic energy dispersion 



parallel to the interfaces is assumed. The field quantizes 
the motion of the electrons in the xy plane, giving rise 
to Landau levels with energies e„ = [n + l/2)Tuo c , where 
n = 0, 1, 2, .. is the Landau index and lo c = eB/m*c is the 
cyclotron frequency. To a good approximation the longi- 
tudinal degrees of freedom are decoupled from the trans- 
verse motion and may be treated independently. The 
probability amplitude 6™" for an electron in a time de- 
pendent state | a > to be at plane j along z and in the 
Landau level n is determined by the equation of motion 



dh na 

in ^- = (e j + e n + Uyjbff)b^ 



dt 



m/3 



(i) 



v(b™ 1 +b™ 1 ~2b™) 1 



where v is the hopping matrix element between electrons 
in nearest neighbor planes, and ej represents the band 
contour and external bias. Here the sum over (m, 0) cov- 
ers all occupied electron states of the system, within the 
various Landau levels m with energies below the Fermi 
level, incident on the DBS from the emitter side. In 
writing Eq. (1) we have adopted a Hartree model for 
the electron-electron interaction, keeping just the intra- 
atomic terms as measured by the effective coupling con- 
stant U. We have neglected inter-Landau level transitions 
since thjey are of little importance at not too low magnetic 
fields,M and averaged over allowed transitions between 
degenerate states within a Landau level, taking advan- 
tage of the localization induced by the gaussian factor 
in the Landau basis. Detailed inclusion of long range ef- 
fects make numerical calculations more difficult, without 
introducing a qualitative difference in the physical de- 
scription of the system. As we will show in what follows, 
the nonlinear term proportional to U appearing in the 
above equation is of key importance in the behavior of 
the system. 

The time dependent Eq.(l) is solved using a half- 
implicit numerical method which is second-order accu- 
rate and unitaryO. Discretizing the time variable, it can 
be written in the form, 
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\ (^(t + *«) + &?£(* + ft)) (2) 
+ ( J- - 1 - 2«(e, + e n + V?(t)))&?°(* + ft) = 
-!(&??!(*) +&??i(t)) 
+ (-l- + l + 2v(e j + e n + Vf(t)})b™(t), 

where V^(i) = C/£ m(3 l&J^WI 2 and ft is the time incre- 
ment. In solving the tridiagonal matrix of Eq.(2), bound- 
ary conditions must be specified at the left (z = —L) and 
right (z — L) edges of the structure. The approach taken 
here assumes that in the presence of a bias, the rwave 
function at time t outside the structure is given bytil 

b na = j e tk a z 3 + Rjne -ik a z^ z . < _ L (3) 

b™=T jn e< z *, Zj >L. (4) 

Here k a and k' a = \J1ra*\t a — ei]/h are the wavenum- 
bers of the incoming and outgoing states, respectively, 
with e a — e n — Av sin (k a a/2) the energy of the incom- 
ing particle. To model the interaction with the particle 
reservoir outside the structure the incident amplitude I 
is assumed to be a constant independent of the coordi- 
nates. The envelope function of the reflected and trans- 
mitted waves, Rj n and Tj n , are allowed to vary with j, 
however. Since far from the barriers these quantities are 
a weak function of the coordinate Zj we restrict ourselves 
to the linear corrections only. This approximation is ap- 
propriate provided the time step ft does not exceed a 
certain limiting value. For the results presented here, a 
maximum value of ft = 3 x 10~ 17 s was found sufficient 
to eliminate spurious reflections at the boundary while 
maintaining numerical stability up to 20 x 10~ 12 s. 

In our numerical procedure the coefficients obtained for 
one bias are used as starting point for the next bias step. 
Once 6™ Q are known, the current at site j is numerically 
obtained from0 

Jj = M^n?T.j k Q n> imvrwi - 6 D>*- ( 5 ) 

where l m = (hc/eB) 1 / 2 is the magnetic length, k n f = 
yj2m*{tf — e n )/h with e/ the Fermi energy, and the sum 
is over magnetic energies e„ < e/. 

We next apply our model to an asymmetric 
GaAs/AlGaAs double barrier structure, with emitter 
and collector barrier thicknesses of 1.12 nm (2 sites) and 
3.36 nm, (6 sites) respectively, and a well thickness of 
11.2 nm (20 sites). The second barrier is made wider than 
the first in order to enhance the trapping of charge in the 
well. For this geometry the first resonance at zero bias 
and magnetic field occurs at 30 meV. The conduction 
band offset is set at 300 meV. The buffer layers are uni- 
formly doped up to 3 nm from either barrier, so as to give 
a neutralizing free carrier concentration of 2 x 10 17 cm~ 3 



at the contacts. In equilibrium and at B = T, the 
Fermi level lies 19.2 meV above the asymptotic conduc- 
tion band edge, so that the zero bias resonance lies well 
above the Fermi seapJhe parameter values in Eq. (1) are 
set at v = -2.16eW3 and U = lOOmeV. The latter was 
chosen phenomenologically so as to fit the experimental 
I-V characteristic for a GaAs devices in the absence of an 
external magnetic field.Q. The sample has 400 sites and 
the normalization of the wave functions is chosen so that 
charge from the electrons filling up to the Fermi .energy 
exactly cancels the positive charge at the contactsE3. We 
solved Eq. (2) using the procedure described above, for 
an energy mesh appropriate to compute the integral in 
Eq.(5). Good convergence was found for a mesh of 100 
points. 

For small values of the magnetic field a stationary 
solution is reached after a short transient. At large 
enough bias two stationary solutions emerge, however, 
reached separately by either increasing or decreasing the 
applied voltage. In the former case the solution sustains 
a charged well while in the latter the well is uncharged 
and the current is very small. The effects mentioned 
have been obseHied experimentally^ and are presently 
well understood.!!!!! 

A completely novel feature starts to develop as the 
magnetic field is increased. At small values of the exter- 
nal bias a stationary solution is rapidly reached. As the 
bias is increased, however, a critical value arises beyond 
which no stationary solution exists, and the system be- 
gins to oscillate in a self-sustained way. After a range of 
voltages over which the oscillations persist a stationary 
solution is reached again. This is illustrated in Fig. 1 
where the I-V characteristic for an applied field of 13T 
is exhibited. The lobe in the figure shows the current 
maxima and minima in the unstable region, its size and 
location depending sensitively on the magnetic field. For 
the case of Fig. 1 its onset is at a rather low bias and 
reaches into the bistable region, while at lower magnetic 
fields the lobe is entirely contained in the low bias range, 
and lies outside the region of bistability. In fact, our 
results show that the well know bistability and the mag- 
netic field induced instability are entirely separate phe- 
nomena. 

In Fig. 2 we show the current at the center of the well 
for three different values of the magnetic field and a fixed 
bias V = 0.27 V. For B = 9 T, the system is seen to 
reach a stationary state after a transient lasting about 
5ps (Fig. 2(a)). For B = 13 T however, the transient is 
followed by an oscillation that is never damped out (Fig. 
2(b)). Although not perfectly periodic, the oscillation 
has two strong Fourier components at frequencies v ~ 0.3 
and 0.8 THz. The reason why the presence of two fre- 
quencies is predominant remains unclear. As shown in 
Fig. 2(c) for B = 17 T, at still higher magnetic fields the 
oscillation becomes irregular. A power spectrum of this 
latter signal shows that the narrow peaks observed below 
1 THz are replaced by a broad low frequency structure 
suggesting that a chaotic regime has been reached. Note 
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that the initial condition from which all three curves were 
obtained is the stable solution at 0.2 V. When reaching 
a bias of 0.27 V from a closer voltage the transient oscil- 
lations are less pronounced, yet their relaxation time is 
of the same order as in the more extreme voltage jump 
exhibited in the figure. Although the amplitude of the 
sustained oscillations depends somewhat on the initial 
conditions, their frequency and structure does not. 

Figure 3 shows the regions where the various structures 
we have described appear in the two parameter space V- 
B. The thick lines delimit the range in which two stable 
solutions may exist for the same applied potential. The 
dotted regions mark areas where self-sustained oscillatory 
solutions appear. Regions with highly irregular oscilla- 
tions suggesting a chaotic behavior appear filled with a 
brick design. We note that while we include in our sim- 
ulations up to several hundred thousands time steps, the 
number of oscillations covered before numerical instabili- 
ties arise are not enough to obtain a fully detailed power 
spectrum at all values of parameters. The search of sta- 
tionary solutions using a standard self-consistent itera- 
tive loop, however, did show unambiguously a transition 
to chaos through successive bifurcations within the re- 
gion where no stationary solutions exist. Although the 
iteration number in this case is not a true time variable, 
the parameter values at which non stationary solutions 
exist and the nature of the latter coincide with those ob- 
tained in the true-time dependent analysis. We used this 
fact in sketching Fig. 3. 

We interpret our findings in the following way. In gen- 
eral, current flows through the system as long as a tun- 
neling resonance lies within the emitter Fermi sea. The 
resonance acts effectively as an energy filter for the trans- 
mitted electrons. We assume that at very low bias it lies 
above the Fermi energy so that no current flows. As 
the bias increases the resonance drops, reaching eventu- 
ally the Fermi energy, thus opening a channel for elec- 
trons to tunnel through the double barrier. A current is 
thus established. What happens next depends sensitively 
on the presence or absence of a magnetic field. In the 
absence of the electron-electron interaction the current, 
roughly speaking, grows linearly with the bias if B=0. 
This behavior is due to the gradual increase in tunneling 
states that satisfy the conservation laws as the voltage is 
raised. At a critical bias the falling resonance reaches the 
bottom of the conduction band, after which the current 
drops abruptly. This overall behavior gives the I-V curve 
its characteristic triangular shape. For finite B however, 
also in the absence of electron-electron interaction, owing 
to the new density of states each time a Landau level en- 
ters the Fermi sea the current rises abruptly, remaining 
essentially constant until the next level comes in. The I- 
V curve acquires then a staircase shape, with a step rise 
made smooth by the resonance line shape. When the 
resonance reaches the bottom of the first Landau level 
in the emitter, the current falls abruptly, as in the B=0 
case. 

We next take into account the electron-electron inter- 



action. As the resonance enters the Fermi sea, current 
begins to flow and charge is trapped in the potential 
well, rising its bottom and pushing the resonance towards 
higher energies. The current drops, some of the accumu- 
lated charge leaks out allowing the current to flow more 
easily once again, and a new cycle begins. When B=0 
this process occurs at very low current, and as we have 
verified numerically, the ensuing oscillations are damped 
out in all cases. At finite field however, as the resonance 
enters the Fermi sea, the abrupt rise of the current to a 
large value triggers oscillations that may be sustained by 
the feedback mechanism if the field is large enough. 

The transition from damped to sustained oscillations 
may be appreciated in Fig. 2. As the magnetic field is 
raised, the degeneracy of the resonance increases thus al- 
lowing more charge to be trapped between the barriers. 
The nonlinear coupling in Eq. (1) thus grows with the 
field, making the latter an effective tunable parameter 
for the amount of nonlinearity in the system. Note also 
that while the dynamic bistability characteristic of DBS 
is related to the drop of the resonance below the bottom 
of the Fermi sea, the effect we are discussing arises as the 
resonance enters the Fermi sea, in support of our numer- 
ical finding that they are entirely independent effects. 

According to the picture drawn above the charge in 
the well lags the current, as exhibited in the inset of Fig. 
2(b). Here the charge at the center of the well (dashed 
line) is plotted together with the current at the same 
point (full line), the former displaced a time r ~ lAps to 
the left with respect to the latter. The relaxation time 
r and the period of the oscillations are determined by 
the tunneling time for the electrons to leak out through 
the barriers and may be adjusted by modifying the bar- 
rier thicknesses. The oscillatory current is always non- 
uniform accross the device, with prominent oscillations 
in the well and the emitter side, while the collector cur- 
rent exhibits a weaker structure. This is due to the 
large width of the collector barrier relative to that in 
the emitter side, a necessary feature of design to have 
enough charge accumulation in the well, and thereby, a 
strong electron-electron interaction. In fact, nonlinear 
effects disappear in our sample if the width of the col- 
lector barrier is reduced from six to just thcee sites, in 
agreement with experimental observations Jiil The very 
high frequency oscillations seen early in the process and 
shown in the inset of Fig. 2(c), are quickly damped out 
and do not contribute to the long term behavior of the 
system. We attribute them to a transient pulse going 
back and forth in the space between the barriers before 
it looses coherence and decays. Its period corresponds 
to electrons at the Fermi velocity being successively re- 
flected by the barriers in the well region. 

There is a wide region of accessible values of bias 
and magnetic field where oscillatory currents are either 
present or absent. This opens up the interesting and 
unique possibility of studying experimentally the tran- 
sition from stationary to oscillatory behavior and possi- 
bly chaos, as the magnetic field is increased in an asym- 
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metric DBS device under bias. When oscillating, the 
system could act as an electromagnetic generator in the 
THz region. Research incorporating the interaction be- 
tween the current and the radiation field is currently un- 
der progress. 

Work supported by FONDECYT 

grants 1950190, 3950026 and 1960417, Fundacion An- 
torchas/Vitae/Andes grant B-11487/4B005, CNPq and 
FINEP. 



FIG. 1. I-V characteristic of the asymmetric DBS at 
B = 13 T. The lobe edges mark the oscillation extremes 
when an unstable solution is present. The dashed line repre- 
sents the lower arm of the bistable regime, reached when the 
bias is being decreased. 

FIG. 2. Current as a function of time at the center of the 
well for V — 0.27 V and an applied magnetic field of (a) 
9T, (b) 13T and (c) 17T. The inset in (b) shows the charge 
and current shifted one relative to the other by 1.4 ps. The 
inset in (c) exhibits a detail of the high frequency transient 
oscillations. 
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FIG. 3. Phase diagram in the B,V parameter space. The 
thick lines delimit the region where bistability occurs. Dot- 
ted areas correspond to oscillatory behavior, while the brick 
design marks the regions where highly irregular solutions are 
present. The filled circles mark the values used in Fig. 2 
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